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Abstract 
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order. The result is found to be consistent both with the expected renormalon 
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1 Introduction 



In recent years, the Numeric implementation of Stochastic Perturbation Theory (NSPT) was 
introduced, which was able to reach unprecedented high orders in pertubative expansions 
in Lattice Gauge Theory (LGT). The historic success of NSPT was the computation of the 
Lattice SU(3) basic plaquette to order /3 _8 |l|]. This led to the possibility of actually veri- 
fying the expected dominance of the leading InfraRed (IR) Renormalon associated to a 
dimension four condensate. The computation was performed on a 8 4 lattice. We are now 
in a position to quote results to a reasonable accuracy for the perturbative expansion of the 
same quantity to order (3~ 10 on both a 8 4 and a 24 4 lattice. 

There are at least three main motivations for such a computation. First of all, it provides a 
clear consistency argument in support of the conclusions of our previous work. The results 
in are actually the only ones from which one can recognize the presence of renormalons by 
direct inspection of the perturbative coefficients of a QCD quantity. Ironically, this is true 
in a scheme (the lattice) which one would never choose for a high order QCD calculation 
in a standard (as opposed to NSPT) approach. We will give strong evidences that a very 
good control has been taken over both the IR renormalon growth and the finite size effects 
on top of that. In particular, finite size effects have always been a matter of concern for the 
computation. They play a fundamental role since they impinge on the IR region which is 
responsible for the renormalon growth. Basically, the higher is the loop, the lower are the 
(IR) momenta which give the main contribution. In view of this one could be worried that 
one is missing important contributions on the 8 4 lattice on which the NSPT perturbative 
computations were first performed on. In |J evidences from the non linear sigma models 
were collected which suggested that renormalon growth was not substantially tamed at the 
eight loop level. We can now strengthen our confidence that the leading behaviour had been 
actually properly singled out in our previous work. On top of that we are now in a position 
to assess the finite size effects, which turn out to be of the order of a few percent at the order 
we got. 

While the eight loop expansion in [|]] was a major success, it actually opened more problems 
than it settled down and this sets the stage for the second motivation for our present com- 
putation. There was a longstanding problem concerning the possibility of determining the 
Gluon Condensate (GC) from LGT by a method which actually relied on the computation 
of high orders in the perturbative expansion of the plaquette. The latter is the lattice rep- 
resentation for the GC and its perturbative expansion is an additive renormalisation which 
is by far the leading contribution at the value of (3 usually taken into account. This per- 
turbative contribution can be seen as the contribution coming from the identity operator 
in the Operator Product Expansion (OPE) for the plaquette. The key point is that in this 
OPE there is no space for an order a 2 term, a gauge invariant operator of dimension two 
missing. Next term in OPE is expected to be an a 4 contribution associated with the genuine 
GC. The old approach to the problem was to subtract the first (perturbative) contribution 
from Monte Carlo measurements of the plaquette in order to single out the a 4 term (the 
GC) [|]]. The signature of such a term can be recognized by performing the subtraction 
at different values of /3 and looking for the scaling dictated by Asymptotic Freedom. The 
actual problem is that the series for the plaquette "not only has to be subtracted, it has 
to be defined in first place"! 2 ]. Due ^° ^ ne P resence °f the IR renormalon, the series has an 

2 This is a citation from Beneke's review 
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ambiguity in different prescriptions for resumming it which is just the same order as the GC, 
i.e. an a 4 contribution. Being aware of this, the most striking point is that by performing 
the subtraction, with any given prescription, one is left with an a 2 contribution, at least 
on the 8 4 lattice the computation was first performed on |J. The problem challenges our 
understanding of potentially any short- distance expansion and so it needs to be cleared up. 
At the moment more than one explanation appear to be possible. The effect could either 
be a pure lattice artifact or a fundamental issue. With respect to the latter hypothesis, 
suggested the effect of Ultra Violet (UV) contributions to the running coupling. Since the 
time of || the phenomenology of non-standard Q~ 2 power effects has actually been reported 
growing. In particular, the work in || establishes a relation between a dual superconductor 
confinement mechanism and "OPE-violating" corrections. In this paper we do not want to 
enter at all the a 2 -affair. Still, the reliability of our previous perturbative computations has 
to be assessed with this question in mind. Given the potential impact of the result in [H , we 
want to be sure that nothing was wrong with respect to the Renormalon contribution extrac- 
tion, leaving for a subsequent publication M a refinement of the whole analysis concerning 
the subtraction procedure. We stress that such a refinement will in particular benefit by the 
study of finite size effects. 

Finally something can be learned from the computation about the NSPT method itself. The 
present result establishes the highest order ever reached by NSPT in LGT. In [|[ a word of 
caution was spent on uncritical application of the method, stressing the issue of a careful 
assessment of statistical errors. Ironically enough, it was pointed out that problems plague 
simple models much more than complex systems like Lattice 577(3). We refer the reader to 
a future publication || for a more organic report on the application of the NSPT method 
to LGT. Still, it is reassuring (although of course admittedly not a proof of reliability of 
the method) to note that in the context of this work there is a fine agreement between the 
results obtained by our stochastic method and those inferred by a (maybe strong) theoretical 
prejudice. 

The paper is organized as follows. In Sect. 2 we recall the basics of IR Renormalon analysis 
applied to a dim = 4 condensate and write the formulas relevant for a finite lattice, which 
are the foundations of our subsequent analysis. In Sect. 3 we present our results for the first 
ten coefficients in the expansion of the plaquette both on a 8 4 and on a 24 4 lattice, showing 
that they are consistent with the Renormalon factorial growth and with the expected finite 
size effects on top of that. Sect. 4 contains our conclusions and perspectives for future work. 
An appendix covers algebraic details. 



2 The IR Renormalon on a finite lattice 



We start by writing first of all the definition of our observable, that is the basic plaquette 
normalised as 

W(M) = 1 - i (Tr U p ) , (1) 

Up being the product of links around a 1 x 1 Wilson Loop. We write an explicit dependence 
on the size (M) of the lattice on which the osbervable is computed. In Perturbation Theory 

W(M) = Y,4 M) p- e (2) 
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Next we move to recall the main points of Renormalon analysis. Notations are slightly 
different from those of |5| and closer to those of . The aim of this section is anyway to 
be self-contained. 

2.1 Basics on the factorial growth of perturbative coefficients 

The expected form for a dimension 4, Renormalisation Group invariant condensate is written 
as 

rQ 2 k 2 Hk 2 

W = J Q ^f(k 2 /A 2 ). (3) 

Q is in our case the UV cutoff fixed by the lattice spacing: Q = 71/ a. Given the above 
dimensional and R.G. arguments, f(k 2 /A 2 ) is a dimensionless function independent on the 
scale Q, for large Q, and can thus be expressed in terms of a running coupling at the scale k 2 . 
One obtains the Renormalon contribution by considering the high frequencies contribution 
to Eq. (W), that is 

rQ 2 k 2 Hk 2 

in which f(k 2 /A 2 ) has been taken proportional to the perturbative running coupling (higher 
powers of the coupling simply result in subleading corrections to the formulas we will get 
this way). We now introduce the variable 

z = z (l - a s (Q 2 ) / a s (k 2 )^j , z Q = — , (5) 

using which together with the two loop form for a s (k 2 ) results in 

W ren = Af ~ dz e~ pz (z - z)- 1 ^ . (6) 
J 

In the last equation we have traded a s (k 2 ) for the (3 coupling one is more familiar with on 
the lattice and introduced a couple of new symbols according to 

4vra s (Q 2 ) = 6/(3 , 7 = 2^, < z < z _ = z (l - a s (Q 2 )/a 5 (rA 2 )) . (7) 

% 

In the equations above bo and b\ are the first and second coefficients of the perturbative (3- 
function (to fix normalisation: bo = ll/(47r) 4 ). zq_ is clearly reminiscent of the IR cutoff rA 2 
imposed in Eq. (f|) to avoid Landau pole once the perturbative coupling had to be plugged 
in. From Eq. (^j) it is now easy to obtain a perturbative expansion 

The factorial growth of the series is the Renormalon contribution, while the O(e~ zo13 ) is again 
coming from the IR cutoff r A 2 in Eq. (^) . Naively one could have not imposed any IR cutoff 
before going from the ill-defined (because of the Landau pole) integral to the series. The 
factorial growth would have anyway asked for the O(e~ zo/3 ) contribution in order to obtain a 
resummation of the series. Actually the reader familiar with the subject has most probably 
well recognized in Eq. @ a Borel representation. 
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2.2 The factorial growth on a finite lattice 



Apart from the explicit reference to the lattice spacing in the UV cutoff we have till now 
adhered to a continuum notation, implicitly assuming some continuum scheme in which for 
example Eq. (|3|) is defined. We now move to rewrite such a representation on a finite lattice. 
One can consider (notice that we now explicitly write the dependence on the number of 
lattice points M) 

rQ 2 k 2 dk 2 

W-(M) = C 2 k -^a s {se). (9) 

The lower limit of integration is the explicit IR cutoff imposed by the finite extent of the 
lattice (Qo(M) = 2tt/ Ma), while the scale s is in charge of matching from a continuum 
to a lattice scheme. In view of that and since rescaling k 2 — > sk 2 amounts to rescaling 
A — > s _1 ' 2 A, it is likely to expect for s a value s = K~ 2 consistent (up to an eventual 
further change of scale) with A cont = KAi att for some continuum scheme. By a change of 
variable in the integral it is easy to trade the rescaling of the argument of a s for a rescale of 
Q 2 — > sQ 2 (and similarly for the lower limit of integration). Remember that in the previous 
section we obtained an expansion in the coupling defined at the scale Q: at one loop level 
a change of scale in the coupling is equivalent to a change of scheme, and this matches the 
present notation to the notation of |], Given Eq. @, one can of course perform the 
same change of variable defined in Eq. (|5|) and obtain a new power expansion. In order not 
to have the main message of the procedure obscured by trivial algebra, we collect the details 
in an appendix. One is left with 

W Icn (M) = P~ £ c? n {M; s, C) (10) 
e=i 

The coefficients c r e en (M; s, C) can be expressed in terms of incomplete F functions (see the 
appendix). It is important to notice their dependences. They depend trivially (i.e. multi- 
plicatively) on the overall constant C which is present in Eq. (0) to account for everything 
which is subleading. They also depend on the scale s which is in charge of the matching to 
the scheme in which Eq. (Q) best describes the expansion (remember that one does not know 
that a priori). Finally, they depend parametrically on the lattice size M. The attitude to 
take with respect to Eq. (|ToD is now the same as in [Jj, [5|, that is: can one determine values 
C = C and s = s by fitting the c£ en (M; s, C) to the high orders c\ (i.e. the coefficients 
actually computed on a given lattice size M)? A positive answer states that the perturba- 
tive expansion of the plaquette is actually described by the asymptotic dominance of the 
leading IR Renormalon. Needless to say, once the fit has been done and values s and C have 
been obtained, the c£ en (M; s, C) contain in a sense all the information needed to test the IR 
Renormalon dominance. Given the data at our hand (coefficients c[ M ^ for £ < £ max = 10 and 
M — 8, 24), one can perform a variety of checks. For example, one can fit C — C and s = s 
from data up to £ < £ max and then compare results at higher loops. As for finite size effects, 
one can fit the relevant parameters on a given lattice size M (remember that the dependence 
on M is a parametric one) and then compare c\ en (M' ^ M\ 8, C) with the actual c[ M>) . Such 
comparisons are what we actually proceed to do in the next section. 

Before moving to the actual results, we think it is worthwhile to recall something about the 
crucial issue of finite size effects with respect to the dependence on both M and s of our 

3 For further details see later and the appendix. 
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formulae. The most straightforward way to obtain the IR Renormalon factorial growth goes 
through the most direct exploitation of Eq. (|4]): one simply insert for a s (k 2 ) the leading 
logs (one loop) formula. From this simple recipe one obtains the T functions which are the 
signature of Renormalons. For a dimension 4 condensate (as the observable at hand) order 
/ turns out to be proportional to (we write the explicit dependence on the IR cut-off) 

fQ' dk 2 k* t _ 1( <3 2 , 

At this point a simple steepest descent argument points out that momenta k* that contribute 
the most to the £th order decrease exponentially with i |1 

k* ~ e-^-D/ 2 (12) 

and that's the reason for being a priori worried for a substantial taming of the factorial 
growth on a finite lattice: starting from a certain order momenta k* will fall outside the 
integration region. Now the point is that our formulae actually ask for the scale s in front 
of the k 2 in the argument of the logarithm which results in 

k* ~ s-V^ 1 )/ 2 . (13) 

Being s << 1 this changes quite a lot the estimate of the dependence on M of relatively high 
orders in Perturbation Theory. We will see that finite size effects are well in accord with our 
formulas and under control at the order we got. 



3 The Lattice SU(3) plaquette to tenth order 

In Table (1) we quote the coefficients of the expansion of the 577(3) plaquette to tenth order 
on 8 4 and 24 4 lattices according to the notation of Eq. (0). The results were obtained on the 
APE-systems the Milan-Parma group were endowed with by INFN. Reaching ten order (up 
to the errors one can inspect from the table) required for the 24 4 lattice about 2000 hours 
on a qhl system (128 nodes, delivering a peak performance of 6.4 Gflops). We stress that 
this system has got an expanded-memory; the 4 MWords on each node allow space enough 
for a fairly big lattice to such an high order (we recall that in NSPT each field is replicated 
a number of time which is twice the order in a s , i.e. the order in g). Results for the 8 4 
lattice came out of about 4000 hours on a q4 system (32 nodes arranged in 4 boards, on 
each of which a replica of the lattice was allocated; this system delivers a peak performance 
of 1.6 Gflops). In considering these data and their statistical significance one should keep 
in mind that the signal to noise ratio is better for the bigger volume. By inspecting the 
table one can make some first comments. First of all, we improved the errors with respect 
to 0. To improve errors on low loops would have been relatively easy, but since we were 
mainly interested in high loops, the simulation times (and therefore the statistical errors) 
were actually driven by the tenth ordeiQ. Having said that, we decided we had the statistic 
significance we needed when finite size effects were disentangled from statistical errors. By 

4 The time spent on a single iteration basically scales as (O 2 — 0)/2, O being the order in For 
more details on technical points concerning the application of NSPT to LGT see M . 
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Table 1: The first ten coefficients in the perturbative expansion of the Lattice SU(3) pla- 
quette. 



n 


c 8 
°n 




c 24 




1 


1.9994(6) 


10° 


2.0000(4) 


10° 


2 


1.2206(16) 


•10° 


1.2208(10) 


• 10° 


3 


2.9523(58) 


•10° 


2.9621(48) 


• 10° 


4 


9.345(27) 


10° 


9.417(29) 


10° 


5 


3.397(14) 


10 1 


3.439(20) 


10 1 


6 


1.346(7) ■ 


10 2 


1.368(9) • 


10 2 


7 


5.653(34) 


10 2 


5.774(42) 


10 2 


8 


2.480(18) 


10 3 


2.545(21) 


10 3 


9 


1.124(10) 


10 4 


1.159(10) 


10 4 


10 


5.227(52) 


10 4 


5.416(57) 


10 4 



inspecting the coefficients one can see that finite size effects between 8 4 and 24 4 are less than 
1% for the first four orders, less than 2% for the first six, less than 3% for the first eight 
and less than 4% up to the highest order we got, i.e. tenth order. A proper infinite volume 
extrapolation is out of reach at this level, since one would need better statistic and of course 
more lattice sizes. Still, we will see that within the IR Renormalon dominance approach it 
is easy to estimate the correct order of magnitude of such an extrapolation. As we will see, 
up to tenth order one can infer a correction of the 24 4 results with respect to infinite volume 
of less than 0.5%. 

Having in mind all the points we have already made in the previous sections, we now approach 
the two basic questions that set the stage for this work, that is: Was the IR Renormalon 
growth correctly singled out in |1|, ||? Are finite size effects under control? We will show 
that there is a positive answer for both questions. 

3.1 The leading Renormalon behaviour 

In order to test the IR Renormalon dominance claimed in [[l], [|, we now proceed to perform 
some consistency checks, just of the type we have already referred to in Section ( |2.2| ). These 
come out of three different fits which we now proceed to describe and which are in a sense 
a summary of a variety of checks we performed on our results. 

• First of all we try to understand to which extent the ninth and tenth order (new) 
results could be inferred from the first eight (already computed) coefficients. In order 
to do that, we fit the formula for c r e en (M; s,C) to the set {<4 , cf A \ cf\ cf 4 ^} (notice 
that while we want to make contact with previous results, nevertheless we now take 
into account both lattice sizes). We will refer to this procedure as Fit 1. As the result 
of this fit (and later of the following fits) we quote s x = y/s instead of s. This quantity 
is simply related to the matching between the A parameters between the continuum 
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10= 



4 5 6 7 

perturbative order I 



Figure 1: The comparison between (red crosses) and c r e en (24; s, C) as inferred from Fitl. 
The latter are marked as circles for the orders that have been taken into account in the fit; 
they are instead marked as dots or diamonds for other points. In particular diamonds can 
be looked at as forecasts for the two highest orders computed. 



and lattice scheme one is quite familiar with 



A 



Si = V s — x ■ 



latt 



A, 



(14) 



emit 



In the previous formula the factor a; is a possible further change of scale (see later). 
We obtain a value for the the scale S\ = 0.0074. In previous works by our group 
the fits were performed without referring to a change of scale in the integral. The 
matching from a continuum to the lattice scheme was obtained by a change of variable 
Pcont = Piatt — r — -J—. In terms of the notation of this work the result of the fit in 

Piatt 

[§] reads Si = 0.0078. The 5% difference sets a rough order of magnitude for what one 



Table 2: Perturbative coefficients c r e en (M; s,C) as inferred from Fitl. 5m are per cent 
deviations from the actual values of c\ . 5oo are the residual per cent deviations of 24 4 
coefficients with respect to infinite volume. 





M 


= 8 


h 


M = 


24 


#24 


<$oo 


ren 
C 6 


1.30 


■10 2 


-3% 


1.33- 


10^ 


-3% 


< 0.1% 


„ren 

c 7 


5.71 


•10 2 


+ 1% 


5.84- 


10 2 


+1% 


-0.1% 


„ren 
c 8 


2.51 


• 10 3 


+ 1% 


2.58- 


10 3 


+2% 


-0.1% 


„ren 
c 9 


1.11 


■ 10 4 


-1% 


1.15 • 


10 4 


-1% 


-0.2% 


„ren 
c 10 


4.90 


■ 10 4 


-7% 


5.13- 


10 4 


-6% 


-0.3% 
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has to live with in this and the other comparisons that we make: one should always 
keep in mind that we are managing notations that differ by higher orders. Notice 
that this agreement is a first good message for our understanding of finite size effects 
based on our formula for c r e en (M; s,C): since we now fit results both on the 8 4 and 
on the 24 4 lattices (while the results in || were based on 8 4 alone) the parametric 
dependence on the lattice size M is supposed to work pretty well (anyway see later for 
more definite statements). With respect to Eq ([14]) Beneke |2| made the point that 
the result Si = 0.0078 in || could for example be interpreted by matching to the MS 
scheme and taking into account a further change of scale (the factor x in Eq (|I4"D). The 
latter amounts to computing the MS coupling at a scale 0.706/a instead of n/a. With 
the second (standard) choice of scale one obtains the standard result A^g = 28.8 Ai att . 
Sl = 0.0078 could in turn be interpreted as 0.706/tt x 1/28.8 = 0.0078, while the 
result si = 0.0074 could be read 0.666/vr x 1/28.8 = 0.0074. One should not take all 
these considerations too seriously, e.g. there is no compelling commitment to the MS 
scheme (once again, remember that one does not know which is the scheme in which our 
theoretical prejudice works the best). Still, it is reassuring that the numbers one gets 
are absolutely sensible as for their order of magnitude. Having made contact with our 
previous results, we now proceed to make contact with what we are mainly interested 
in, that is the new (ninth and tenth) orders. One can get a first glance at this by looking 
at Fig. (1) in which we compare the computed cf^ with the c£ en (24; 0.0074 2 , C^) (we 
will not quote the values obtained for the overall constant C since they are not too 
enlightening). Notice that different symbols refer to different orders as far as the 
c£ en (24; 0.0074 2 , C^) are concerned: we want to emphasize the difference between the 
orders that are actually included in the fit and those that are looked at as forecasts to 
assess how well the asymptotic behaviour has been singled out. A very good agreement 
is manifest, which is actually magnified by the logarithmic scale. A careful inspection 
of the numbers themselves can be got from Table (2) in which we quote the values of 
the last five coefficients as inferred from Fit 1 on both lattice sizes. A first point to be 
made has to do with finite size effects. As we have already mentioned, a fit to data on 
both sizes is absolutely sensible and the deviations of the inferred coefficients from the 
actual ones are much the same on 8 4 and on 24 4 . Notice also that the residual finite 
size effects between 24 4 and infinite volume turn out to be negligible. As for the crucial 
issue of how well Fit 1 can single out the asymptotic behaviour we notice that the 
ninth coefficient is forecast with a precision of 1%, while a 6-7% deviation is there for 
the tenth order. Notice also that, while orders seven and eight are kept into account 
in the fit, the agreement for orders less than seven results in a sensible shape for an 
asymptotic behaviour (by the way, the results are pretty stable to the inclusion of order 
six in the fit). This is the right time to make a point that specifies a statement that we 
have already made, i.e. that in inspecting the forecasts of our fits one should keep in 
mind that we are managing expressions that are there up to higher orders. This is of 
course not the end of the story. An equally important indetermination is coming from 
the fact that only a posteriori one can understand how asymptotic the expansion was 
at the highest order which has been taken into account in the fit. Roughly speaking, 
by fitting only the expansion for one quantity there is no obvious way to prevent the 
fit from pretending to have already reached the actual asymptotic regime. 

• With respect to the last point we made there is a nice way to improve. We exploit it 
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in what we call Fit 2. This time we make use of the fact that in [T|] results for the first 
eight orders were computed not only for the lxl Wilson loop (the basic plaquette), 
but also for the 2x2. 



°- 2 
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perturbative order I 



9 10 11 



Figure 2: The comparison between cp (red crosses) and c r e en (8; s,C) as inferred from Fit2. 
Same notations as in Fig. (1). 



Now, the IR Renormalon dominance is supposed to force a universal asymptotic be- 
haviour. That's why we now fit the formula for c r e en (M; s, C^) to the set {c^\ 



(8,2x2) (8,2x2) (8,2x2) 



Cg ' cf Cg ' '}. The cf' 2x2 ^ are taken from [JTJ] , while of course this time a cou- 
ple of different overall constant C^ l > have to be fitted, corresponding to i = 1 X 1 and 
i = 2 x 2. Notice that this time, as an extra test of the control on finite size effects, 
we make use only of 8 4 results for the basic plaquette. The value fitted for the scale 
is now Si = 0.0065. While we are still stopping at eight order in what we take into 
account for the fit, the curvature in the growth of the 2x2 coefficients is opposite to 
that of 1 x 1 and this corrects a bit with respect to Fit 1. As for the 12% change in the 
scale, much the same that has already been told holds: while there is in a sense extra 



Table 3: Perturbative coefficients c r e en (M; s,C) as inferred from Fit2. Same notations as in 
Tab. (2). 





M 


= 8 


Ss 


M = 


24 


#24 


5oo 


„ren 
C 6 


1.23 


•10 2 


-10% 


1.25 ■ 


10 2 


-10% 


< 0.1% 


„ren 
c 7 


5.51 


•10 2 


-2% 


5.64- 


10 2 


-2% 


-0.1% 


c 8 


2.48 


• 10 3 


+0.2% 


2.55- 


10 3 


+0.3% 


-0.1% 


„ren 
c 9 


1.12 


• 10 4 


-0.2% 


1.16 • 


10 4 


+0.2% 


-0.2% 


„ren 
c 10 


5.09 


• 10 4 


-3% 


5.31 ■ 


10 4 


-2% 


-0.3% 
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Figure 3: The comparison between (red crosses) and c r e en (24; s, C) as inferred from Fit3. 
Same notations as in Fig. (1). Notice that this time there are no diamonds. 



information on what has to be understood as asymptotic, it is nevertheless reassuring 
that the change in the coefficients is quite smooth and absolutely sensible as we are 
managing asymptotic behaviours. Not surprisingly, one can inspect from Fig. (2) (this 
time we plot the comparison for coefficients on the 8 4 lattice) and (even better) from 
Table (3) that we have fairly improved the asymptotic behaviour. In particular, devia- 
tions for the forecasts on order ten are reduced to 1-2%, while other nice features like 
the consistency of finite size effects are still there. Also the residual finite size effects 
dependency turns out to be just the same (negligible) order. 



One can of course devise a fit in which all the information at hand is taken into account. 
This is what we proceed to do in what we refer to as Fit 3. This time we fit the formula 
for s CW) to the set {cf , c?\ c?\ cf\ c£\ cf> , cf \ c™, cf> , eff , cf^\ 

c (8,2x2)^ c (8,2x2)-^ are a g a j n changing the players on the ground by including a 

wide range of orders on both lattice sizes, together with information from 2x2 Wilson 



Table 4: Perturbative coefficients c r e en (M; s, C) as inferred from Fit3. Same notations as in 
Tab. (2). 





M 


= 8 


5 8 


M = 


24 


#24 


<$oo 


ren 
C 6 


1.14 


■ 10 2 


-17% 


1.16 ■ 


10 2 


-17% 


< 0.1% 


„ren 

c 7 


5.26 


•10 2 


-7% 


5.37- 


10 2 


-7% 


-0.1% 


„ren 
c 8 


2.43 


■ 10 3 


-2% 


2.49 ■ 


10 3 


-2% 


-0.1% 


„ren 
c 9 


1.12 


■ 10 4 


+0.2% 


1.16 • 


10 4 


+0.4% 


-0.2% 


„ren 
c 10 


5.21 


■ 10 4 


< 0.1% 


5.43- 


10 4 


+0.5% 


-0.3% 
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Loop as well. The resulting scale is now si = 0.0056, with a change with respect to Fit 
2 which is again the same order of magnitude as that obtained in going from Fit 1 to 
Fit 2. The result is remarkably stable with respect to keeping into account only the 
last three orders for the basic plaquette (as for the 2x2 Wilson Loop one can at the 
same time keep into account only the last order available, that is order eight). Results 
of the fit are plotted in Fig. (3) for the 24 4 lattice and summarized in Table (4). Of 
course, this time it does not make any sense to compare with Fit 1 and Fit 2 as for 
the accuracy with which high orders are described. 



In the end, what one can state is that going from Fit 1 to Fit 3 one is not going to jeopardize 
the overall picture: the transition in describing the growth of the coefficients is quite smooth, 
consistent with the asymptotic behaviour being better and better described. Again, one is 
validating the finite size effects as embedded in our formulae. 



3.2 Finite size effects 

From all the previous arguments the impact of finite size effects has already been widely 
discussed. Still, we think it is worthwhile to present the tests of finite size effects also in a 
graphical format. In Fig. (4) we depict something that has to do with finite size effects as re- 
sulting from the results of Fit 2. To be definite, we plot the ratios /c£ en (8; 0.0065 2 , C^ 1 * 1 )) 



1.12 



S_1.04 - 



1.02 - 



perlurbative order I 



Figure 4: The ratios cf^ /c r e en (8; s, C) (red triangles pointing up and red diamonds) and 
/c r e en (2A; s, C) (blue diamonds) as inferred from the results of Fit 2. 



and cf >/c\ e n (24; 0.0065 2 , C^). One can directly see the conclusion that has already been 
drawn in the main discussion of Fit 2. While only the 8 4 lattice has been taken into account 
in the fit, the points for both lattices fall on top of each other, i.e. from a fit to 8 4 and the 
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dependence on M embedded in our formulae we describe basically with the same accuracy 
the coefficients of the 24 4 lattice as well. Something similar is shown in Fig. (5) which is 
the same plot for the equivalent of Fit 3 in which only the 8 4 lattice has been taken into 
account. 



1.16 - 



~ 1.1 - 
O 

0) 

3 

g~ 1 08 - 
~°"1.06- 



1.02 - 



7 8 9 

perturbative order I 



Figure 5: Same quantities as in Fig. (4) as obtained by the results of the equivalent of Fit 
3 in which only 8 4 lattice has been taken into account. Quantities for 8 4 are plotted as red 
triangles pointing up while quantities for 24 4 are plotted as blue diamonds. 



4 Conclusions and perspectives 



We computed the perturbative expansion of the basic SU(3) plaquette on both a 8 4 and 
a 24 4 lattices to /3 -10 order. We are now in a position to strengthen the claim in 0,0, 
i.e. the coefficients growth is consistent with both the leading IR Renormalon dominance 
and with finite size effects on top of that. The latter does not exceed 4% at the order we 
got and the residual finite size effects on the 24 4 results are no more than half a percent, 
again at tenth order. Having verified that our previous analysis in || had neither failed in 
the extraction of the leading asymptotic behaviour nor underestimated finite size effects on 
the perturbative coefficients, these results make even more important to refine the analysis 
contained in the same paper as for the extraction of a quadratic contribution to the lattice 
representative for the Gluon Condensate. This refinement will now benefit by a control on 
finite size effects. This is what we plan to do in a near future , maybe going even through 
a further refinement on errors in the perturbative coefficients. 



12 



Acknowledgments 



The authors are grateful to G. Burgio for many interesting discussions and to E. Onofri and 
G. Marchesini for their constant interest and encouragement. F. D.R. acknowledges support 
from both Italian MURST under contract 9702213582 and from I.N.F.N. under i.s. PR11. 
L. S. acknowledges support from EU Grant EBR-FMRX-CT97-0122. 



Appendix 

We now sketch the steps to go from Eq. (Q) to Eq. ([T0| ) . As we have already said in Sec. (|2.2j ) 
one way to treat the dependence on s is to rescale the integration variable (define t 2 = sk 2 ). 
In this way one ends up with the rescaling Q 2 — > sQ 2 . This in turns means (3 — > (3{sQ 2 ) 
which makes contact with the formulation of |l], |5| (again, as already said): (3 — > (3 — r (at this 
level one obtains a one loop formula). We chose another way to proceed which manipulates 
the integrand in another way. By going through the change of variable of Eq. @ one obtains 
(up to overall constants) 

Jo 10 J l + 66 ^ln S ' 

Zi r (M) being the value for z pertaining to the IR cut-off Qo(M). The last factor is simply 
a ^^2y By expanding the latter in a geometric series (we are aiming at an expansion in (3~ l ) 
and by performing a change of variable we end up with 



, ^ f f3z tr (M)-l3z 

/ gT e -^(-l)- 1 -T^(6fo ^lns) n / e- z z- l -^ n dz 

J-/3z 



Note that the effect of the scale s propagates to high orders much the same way the leading 
logs expansion for a s is responsable for the Renormalon growth, i.e. via a geometric series. 
The last input is now the asymptotic expansion for incomplete T functions (which is useful 
once one splits the previous integral in terms of a sum of incomplete T functions) 

°° dte-H a - 1 « z a - x e- z £ Y ^ X J^\ - z )- k 

fe>o ^(1 a) 

Once this expansion (valid for \z\ — > oo in | argz| < 37r/2) is plugged in, one proceeds to the 
final power expansion in (which is easy to manage for example in Mathematica©). 
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